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We investigate neutron-antineutron oscillations in the Randall-Sundrum warped extra dimen- 
sional scenario. The four dimensional effective strengths of the relevant operators that induce the 
oscillations are calculated up to an arbitrary coupling along with their corresponding enhancements 
due to QCD 1-loop running effects. We find that the AB = 2 operators can be geometrically 
suppressed without fine tuning to within current experimental limits with a warped down four 
dimensional mass scale which can be as low as a fraction of a TeV. 



I. INTRODUCTION 

The standard model (SM) of particle physics has, for 
30 years, enjoyed unequaled success in describing the re- 
sults of particle physics experiments. It is, however, not 
an entirely satisfactory theory due to the fact that it has, 
to date, left many unanswered fundamental questions. In 
particular, it provides no explanation for the many differ- 
ent hierarchies which have been built into it. The most 
famous of these being the electroweak-Planck hierarchy 
problem in which, due to the ultra-violet sensitivity of the 
Higgs mass, the massive separation between the Planck 
scale and the electroweak scale is considered to be unnat- 
ural. One particularly appealing solution to this problem 
is the Randall-Sundrum (RS) model [IH3] ■ Within the 
context of this model the large hierarchy arises due to 
the warping of a compactificd Anti-deSitter (AdS) extra 
dimensional geometry. This warping manifests itself as a 
warping factor which exponentially suppresses the mass 
scales within the theory, creating an effective hierarchy. 
Another appealing feature of the RS model is its ability 
to explain the SM fermion mass hierarchies with the same 
mechanism which explains the electroweak-Planck hier- 
archy [SHZ] • By promoting all SM fermions to bulk fields 
the fermion mass hierarchies are explained in terms of 
the fermion geography within the warped extra dimen- 
sional space. In such a scenario, the five dimensional 
(5D) fermion fields are Dirac fields whose wave function 
localization in the extra dimension is completely char- 
acterized by a single 0(1) c parameter. By using a Z 2 
orbifold projection or equivalently by choosing appropri- 
ate boundary conditions on the UV and IR branes one 
can project out the chiral zero modes. The SM fermions 
are identified with these chiral zero modes of the bulk 
fermions and they have exponential wavefunction profiles 
in the extra dimension. The effective Yukawa couplings 
depend heavily on the wave function overlap of the corre- 
sponding fermion wavefunctions with the Higgs, which is 
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situated on the TeV brane in the extra dimension. Heavy 
fermions are localized near the IR brane and thereby have 
a large overlap with the Higgs field, while light fermions 
are localized closer to the UV brane. In this particu- 
lar class of RS flavour models the SM gauge symmetry is 
promoted to a bulk symmetry. The 4D Yukawa couplings 
can then all be taken to be 0(1) while the c parameters 
can all be determined by fitting the fermion masses and 
their mixing parameters [5]. 

One can then ask about the nature of higher mass di- 
mension operators within the context of the RS model, 
such as those corresponding to proton decay and neutron- 
antineutron (n-n) oscillation I9HTT] . If no extra sym- 
metry forbids these operators they will be suppressed by 
some mass scale which is close to or exceeds the UV com- 
pletion scale of the RS model. If one simply takes this 
to be the Planck scale then this would be sufficient to 
satisfy the experimental constraints however the exact 
same warping mechanism which reduces the Planck scale 
to the electroweak scale acts to reduce this mass scale 
suppression as well. 

It is well known that proton decay is a problem for the 
RS model [5]. In order to properly suppress the relevant 
operators for proton decay it is necessary to maintain 
large separations between the quarks and leptons in the 
extra dimension however successful mass configurations 
for these fields do not allow for such large separations. 
As opposed to accepting unnaturally small dimension- 
less couplings for these operators it is thought that there 
exists an extra symmetry which will forbid these opera- 
tors entirely however the exact nature of this symmetry is 
as yet unknown. The simplest solution is to introduce a 
Ux (1) symmetry where X could denote the total baryon 
number (B), lepton number (L), or their difference (B-L) 
which is currently understood to be only an accidental 
symmetry of the SM gauge group. Discrete symmetries 
of the Zm type have also been suggested. 

It is thought that n-n oscillations could present yet 
another problem for models of this type due to the fact 
that the corresponding operators contain only quarks of 
similar mass scales and therefore similar localizations 
within the extra dimension 10 . In the current work 



we analyze the effective strength of the six quark 
operators which induce n-n oscillations in the warped 
RS model assuming that there is no symmetry which a 
priori forbids them. For example, the introduction of 
a or Z3 symmetry would have no effect on the 

operators which induce n-n oscillations but does forbid 
the operators which induce proton decay. Discrete 
symmetries have also been used to study Dirac neutrinos 
in warped models [15] . Previous investigations of the 
effective strength of n-n oscillation operators within 
the context of the 6D Arkani-Hamed-Dimopoulos-Dvali 
(ADD) have yielded a lower bound on the mass scale 
suppression in the observable range Mx ft (45 — 100) 
TeV [13]. 

The paper is organized as follows: Section [IT] briefly 
reviews the treatment of fermions on the 5D AdSs back- 
ground and the problems with proton decay in the RS 
model. Section |III| introduces the relevant operators 
which induce the n-n oscillations and presents the corre- 
sponding 4D effective Wilson coefficients. In section IV 
we present the calculated enhancement of the strength 
of these coefficients due to SM QCD 1-loop renormaliza- 
tion group (RG) running effects which are expected to 
be larger than other gauge interactions and in section [V] 
we discuss our results. 



where E A = diag(e cr , e a , e a , e CT , l/r c ) is the inverse 
fiinfbein, j a = (j^jij 5 ) and m sgn(</>) is the mass. The 
mass must be dependent on the position within the extra 
dimension in order for the mass term to remain invariant 
under the action of the Z2 orbifold. Due to the diago- 
nal nature of the metric the spin connection term does 
not contribute to the action |14j . Any gauge interactions 
are included by simply replacing 8a with the relevant 
covariant derivative. 

The normalized KK mode expansion of the chiral 
fermion fields is chosen to be 
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where the eigenf unctions are orthonormal such that 



(4) 



All SM fields are associated with the zero modes of 
the expansion for which the normalized eigenfunctions, 
determined from solving the zero mode field equations, 
are given by 



II. FORMALISM 

This section serves to define our conventions and no- 
tation. The 5D space is mapped by coordinates x A = 
(x^ , <fi) where the fifth dimension is compactified with size 
r c and <f> G [— tt, tt]. In order to embed a 4D Minkowski 
spacetime within a slice of 5D anti-deSitter space (AdSs) 
with curvature k the points (x^,(f>) and {x^,—^) are 
identified. This creates an S 1 /Z2 orbifold with fixed 
points (2^,0) and (x^,ir). The metric components of 
the warped non-factorizable geometry are given by the 
line element [T] 



ds 2 = G AB dx A dx B = e-^^n^dx^dx" - r 2 c d(\> 2 (1) 

where x^ are the coordinates on the four dimensional 
hypersurfaces of constant <j> with the Minkowski metric 
?7 M „ and <j{4>) — kr c \<p\. Two three branes, called the 
ultraviolet (UV) and the infrared (IR) branes, are placed 
at the orbifold fixed points (f> — and <fi — n respectively. 
The parameters k and l/r c are assumed to be on the 
order of the Planck scale while the product kr c is chosen 
to be ~ 12 to solve the hierarchy problem. 

Working in the low energy effective field theory we can 
write down the free field action for a massive fermion 
ty(x,4>) in the RS background as 
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S = Id xd(j)\ 



kr c (l/2±c LiR ) {lj2±CLR)a 



e 2fcr c 7r(l/2±c iiR ) _ -y 



(5) 



where the normalization factor insures a canonically 
normalized kinetic term in the four dimensional effective 
theory. It is exactly the cl,r = m/k parameters which 
control the localization of the fermion wavefunctions in 
the extra dimension, i.e; cr < 1/2 (cr > 1/2) corre- 
sponds to a closer proximity to the UV (IR) brane for the 
right handed zero mode while cl < —1/2 (cl > —1/2) 
corresponds to a closer proximity to the IR (UV) brane 
for the left handed zero mode. With the Higgs localized 
at the IR brane the effective SUl(2) x Uy(l) invariant 
fermion-Higgs Yukawa interaction with order one cou- 
plings, i/ij, in 5D is given by 



(2) 



d A xd6VG8(6-7r) ^V iL (x, <b)$(x)^ m (x, <f>) + h.c. 
\kr c 

(6) 

where ^ l and are the fermion SUl (2) doublet and 
singlet respectively and $ is the Higgs SUl(2) doublet. 
Integrating out the extra dimensional dependence yields 
the effective mass matrix 



Ma = y^Af^e^+^n) (7) 

where v w = ve~ kTc7T = 247 GeV and we have defined 
the normalization factor 
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equation involving the simple 2x2 Hamiltonian 
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for notational brevity. The nature of this effective mass 
matrix is such that the heavy fermions must be inter- 
preted as being localized near the IR brane and therefore 
have a large overlap with the Higgs held while the light 
fermions must be localized near the UV brane yielding 
a small overlap with the Higgs. Another consequence is 
that only 0(1) differences in the c parameters are then 
needed to generate the observed large effective fermion 
mass hierarchies without any hne tuning of the dimen- 
sionless Yukawa couplings. In particular, using only 0(1) 
differences between the nine c parameters in the quark 
sector, it is possible to reproduce the entire set of quark 
masses and mixing angles. To this extent, there have 
been a number of parameter sets put forward which fit 
the observed data p iT5lH6] . 

Since the RS model is itself considered to be an effec- 
tive theory one can introduce higher dimensional opera- 
tors in the same vein as in other extensions of the SM. 
Some of the more well known examples of these types of 
operators are those responsible for proton decay such as 



d xc 



(giQQQL + g 2 U c U c D c E c ) (9) 



where the fields Q, L, U, D, and E are the bulk ver- 
sions of the corresponding SM fields. Since these oper- 
ators are understood as arising from physics above the 
RS UV cutoff one can conservatively take the mass scale 
suppression M to be on the order of the Planck mass. 
Integrating out the 5D degrees of freedom reveals the ef- 
fective strength with which the zero modes of the above 
fields will induce proton decay. Although the effective 
strength of these operators receives a suppression from 
the resulting exponential overlaps between the various 
fields the Planck scale suppression is warped down and 
replaced by Me~ kr7r . The end result of these two com- 
peting effects is that there is not enough suppression from 
the resulting wavefunction overlap to prevent proton de- 
cay within the current limits and we must therefore con- 
cede either the fine tuning of the dimensionless couplings 
or the introduction of a convenient symmetry (such as to- 
tal lepton number or the above mentioned Z 3 ) which will 
forbid the operator entirely [SJ |S] . Neither of the above 
mentioned symmetries forbid the operators which induce 
n-n oscillations and the question of whether or not the 
RS model can inherently provide the needed suppression 
for these transitions from geometry alone is the subject 
of the present work. 
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where 8m — (n\'H e ff\n} parameterizes the underlying 
physics describing the oscillation [TOl EJ [T7] . The prob- 
ability of finding an antineutron after some time t is then 
given by 



\(n\n(t)}f 



48m 2 



AE 2 + 48m 2 



\/AE 2 + A8mH ) (11) 



where AE = E n — En is the energy splitting. Ex- 
perimental limits from reactors and matter instability 
have produced limits on the off-diagonal components 
of the Hamiltonian of 8m < .75xl0 -32 GcV -1 and 
8m < .6xlCT 32 GeV" 1 respectively [T8H20]. 

The effective Hamiltonian is given by T-L e f / = 
J2i giOi(x, <f>) where each effective operator, Oi(x, <j>), is a 
SU C (3) x SUl(2) x Uy(l) gauge invariant six-quark op- 
erator. Any general AB = 2 operator which contributes 
to n-n oscillations and is constructed from non-scalar 
Lorentz invariant quark couplings can be converted to 
an equivalent operator constructed strictly from scalar 
Lorentz invariant quark couplings via Fierz transforma- 
tions [21j [22] . There are four linearly independent oper- 
ators of this type which are given by [13] 
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03 = (Ql aT CQf )(u R T Cd R ){d R T CcT R )e ij T^ SXT (14) 



04 = {QT T CQ^)(Q^ T CQ l l)(d R T Cd R )e l3 e kl T^ SXT 

(15) 

The round brackets are meant to imply the contrac- 
tion of spinor indices, C is the charge conjugation oper- 
ator, and Greek and Latin indices represent SU C (3) and 
SUl(2) degrees of freedom respectively. The color ten- 
sors contract the SU C (3) indices into color singlet combi- 
nations in two different ways given by [23] 



III. n-n OSCILLATIONS 

The time evolution of an initially slow moving beam 
of neutrons is described by the following Schroedinger 



(16) 



(17) 
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where the first tensor is symmetric about the inter- 
changes (a,0), (7,5), (A,t), (a/3, 7<5), (a/3, At), (■yS, At) 
while the second is anti-symmetric about the inter- 
changes [a, 0\ and [7, 6] and symmetric about (A, r) and 
(a/3, 7<5). These operators can all be easily generalized to 
5D by replacing the fermion fields with the corresponding 
bulk fields. In 5D the coefficients associated with each 
operator have mass dimension -7 so we can rewrite them 
G l 

as Oj = — -pr where the CVs are the dimensionless Wil- 
y M\ 

son coefficients and Mx is the 5D mass scale at which 
a detailed knowledge of the underlying physics respon- 
sible for the generation of the operators becomes indis- 
pensable. In order to obtain the effective dimensionless 
couplings and mass scale in the 4D theory we integrate 
out the extra dimensional dependence as 



The equality of and is due to the fact that 
the two operators 0\ and O2 share the same overall quark 
content and differ only in the way in which the spinor and 
color indices are contracted. Furthermore, as the UV 
complete theory is not known the dimensionless Wilson 
coefficients Ci are also unknown. With no loss of gener- 
ality we then set all the dimensionless Wilson coefficients 
to unity. 

As previously mentioned, there have been a number 
of numerical fits made to the existing observational data 
taking into account not only the quark masses but also 
the CKM mixing angles as well. In table I we reproduce 
a list of three different representative configurations of c 
parameters for the first generation of the left and right 
handed quarks which fit the data [5] 




G\Oi(x, 



eff 



Oi(x) (18) 



where C e J^ is the 4D effective Wilson coefficient 
associated with the i th operator and M 4 o = Mxe~ kr ° v 
is the warped down 4D mass scale which is independent 
of the localization of the particle content in the extra 
dimension. The physics that generates these operators 
is largely unknown and highly model dependent. If 
they arise due to a higher dimensional grand unified 
theory then the mass scale Mx is expected to be of 
order k, the AdSs curvature scale, although if they 
arise due to physics beyond the RS UV cutoff then the 
mass scale can be of order l/r c . In our effective theory 
approach it is appropriate to parameterize the mass 
scale as Mx — pk with p taken to be a free parameter 
of the theory to be determined by experiment. The 
warped down 4D effective mass scale is then written 
as M±d = pke~ kr "' n . This is a convenient choice of 
parametrization since studies of precision electroweak 
measurements and flavour changing neutral currents 
imply that the lowest allowed value of the warped down 
curvature scale is ke~ kTclT = 1.65 TeV [24] . 

Writing the zero modes of the SM quark bulk fields as 

q LiR (x, cj>) = Vk Mlm q L>R {x) e ^ 2±c ^ (19) 
yields the following set of effective Wilson coefficients 
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, ff _ C A Aq l Nj R e krM3+ic QL -2c dR ) 
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TABLE I. Numerical Fits of the Quark Masses with CKM 
Mixing Angles 





Configurations 


c parameters 


I 


II 


III 


c Ql 


-0.634 


-0.629 


-0.627 


Cu R 


0.664 


0.662 


0.518 




0.641 


0.58 


0.576 



where, due to the SUl(2) gauge symmetry, c Ul and 
c dh are equal and we have denoted them both simply as 
cq l . The matrix elements of the 4D effective six quark 
operators (n\Oi(x)\n) have been calculated within the 
context of the MIT bag model in reference [21]. Aver- 
aging the results of the various fits we obtain the val- 
ues (n\Oi(x)\n) = -5.945 x 1(T 5 GeV 6 , (n\0 2 (x)\n) = 
1.485 x 1CT 5 GeV 6 , (n\0 3 (x)\n) = -2.95 x 1CT 5 GeV 6 , 
and (n\0 4 (x)\n) = 2.22 x 10" 5 GeV 6 . 



From Eq.(|18|) the matrix element (n[H e f f\n) will in- 



volve the Wilson coefficients C^- which are determined 
at the scale of M 4 jj whereas n-n oscillations take place 
at the neutron mass scale 1 GeV. This requires us to use 
the relevant RG equations to run the Wilson coefficients 
down to the oscillation scale. The largest contribution to 
the running of the Wilson coefficients will come from the 
QCD sector. This calculation of this contribution is the 
subject of the next section. 



IV. QCD RUNNING OF C\ 



ff 



It is well known that only the SU C (3) coupling runs 
significantly within the range between the TeV scale and 
the neutron mass scale. We therefore restrict ourselves 
to calculating the SU C (3) renormalization effects only. 

g 2 

Working to first order in a s = — we find a total of 15 di- 

47T 

agrams for each operator that must be computed not in- 
cluding the wavefunction renormalization diagrams. The 
first few generic diagrams are depicted in FIG. 1. 
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Permutations 



FIG. 1. Feynman diagrams contributing to the anomaious 
dimension of the effective operators Oi(x) due to SM gluon 
exchange. 



where 7 0i = -Z Q . ^—q^ 



is the anomalous dimen- 



sion of the operator Oi. Through direct calculation the 
counter terms for 0\ and O2 are found to be equal while 
the same is found to be true for C 3 and O4. Using the 
known quark wave function renormalization |25) 



Z a = 1 -In 

q 3tt 



(26) 



the independent operator renormalizations were deter- 
mined to be 



The renormalized effective Lagrangian is 



Za,.. 



IVI 4D i L 

(23) 

The renormalized operators Om. are expressed in terms 
of the unrenormalized operators Oi as 



Oi 



Zrn.O&i 



(24) 



where Ooi are the bare operators which are indepen- 
dent of the renormalization scale fi. Since the renormal- 
ized operator is dependent on p via Zq 1 the 4D effective 
Wilson coefficients must carry a compensating p depen- 
dence to ensure that C e f / is independent of renormaliza- 
tion scale. Since the SU C (3) 1-loop running effects do not 
induce any operator mixing this implies that the effective 
Wilson coefficients each obey a simple RG equation given 
by 



ftr< e ff 
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where A is the RS UV cutoff. These operator renormal- 
izations yield the anomalous dimensions jq 1 2 = and 

4a s (M 4D ) . ' . 

7o 3 4 = where the running coupling is evalu- 

7T 

ated at the intermediate mass scale. Integrating the RG 
equations ( 25 ) down to the the neutron mass scale leads 



to the following scaling behaviour for the effective Wilson 
coefficients 
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(29) 



where m c , to;,, and mt are the masses of the charm, 
bottom, and top quark respectively. The full matrix cle- 
ment, evaluated at the neutron mass scale, which param- 
eterizes n-n oscillations is then given by 



(n\Heff\ 



p 7 (ke- kr - 7T ) 5 



^(GeV 2 )! {n\O x \n) + (n\0 2 \n)) + ( Cf f (GeV 2 )(n\0 3 \n) + C e J f (GeV 2 )(n\0 4 \n) 



~a s {GeV 2 )~ 


8/9 


" Oi s (to 2 ) ~ 


24/25 


~a s {m 2 b )~ 


24/23 


a s (m 2 ) 
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a s (m 2 )_ 




a s (pke~ kr ^) 2 ) _ 





(30) 



where we have factored out the 1 / p 2 dependence from 
all of the effective Wilson coefficients. This leads to a 
simple overall p dependence of the form 



\Sm\ = \(n\H eff \n)\ = - (A + Blnp) (31) 
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where A and B vary depending on which configuration 
is used for the c parameters in the quark sector. In FIG. 
2 1 8m | is plotted for each of the three configurations con- 
sistent with numerical fits to the quark masses and CKM 
mixing angles in TABLE I. 




FIG. 2. 8m is plotted for each of the three configurations 
listed in TABLE I. The horizontal line represents the experi- 
mental limit \Sm exp \. 

For completeness we have included both the curves 
with and without the QCD 1-loop running effects. For 
configuration III the effective Wilson coefficients are such 
that dominates by 2 and 5 orders of magnitude over 
C3 and C4 respectively. The f-loop running effects 
are therefore negligible and the two curves are essentially 
indistinguishable on the present scale. An upper limit on 
p for all three configurations can also be obtained by 
requiring that the value of \8m\ be less than the experi- 
mental limit \8m exp \ = .55 x 10~ 32 [T3"I H"5H2"D] . Imposing 
this constraint for each of the three configurations leads 
to the following bounds: pi > 0.240591, p u ^ 0.568982, 
and pin > 1.96185. We can easily turn these bounds 
on p into bounds on the 4D warped down effective mass 
scale M^n which are then given by M\ D > 0.4 TeV, 
Mlb > 0.94 TeV, and M|j/" > 3.24 TeV. This implies 
that only a relatively small warped down 4D mass sup- 
pression is actually needed to satisfy the currently ob- 
served experimental limits. 



CONCLUSIONS 



We have investigated the effective strength of the lin- 
early independent set of six quark operators which induce 
neutron-antineutron oscillations within the context of the 
warped Randall-Sundrum model. The overall strength of 
the relevant operators arose from the combination of the 
resultant wavefunction overlap of the six quark fields in 
the extra dimensional bulk, the 4D effective warped down 
mass scale suppression, and, to a lesser extent, QCD 1- 
loop running effects. The 4D effective warped down mass 
scale suppression was parameterized by a dimensionless 
factor in order to determine the extent of any extra sup- 
pression needed beyond the minimum allowed by flavour 
changing neutral current constraints. It was determined, 
for the quark c parameter configurations listed, that the 
constraints on the dimensionless factor are such that the 
effective warped down mass scale never has to be greater 
than 0(1) TeV and, in two of the three configurations, is 
only required to be a fraction of a TeV even with enhance- 
ments due to 1-loop running effects. The enhancements 
due to the QCD running were included but were deter- 
mined to not have an overtly large effect on the strength 
of the operators as the contributions from the warped 
geometry far outweighed any running effects. The resul- 
tant wavefunction overlap of the six quark fields in the 
bulk play the most significant role in the suppression of 
the effective operators. The resultant overlap is sensi- 
tively controlled by the c parameters of the quark fields 
which one determines by fitting the quark masses and 
the CKM parameters. The reason that these effective 
operators receive greater suppression than their proton 
decay counterparts stems from the simple fact that the n- 
n transition operators contain more fermion fields which 
leads to more negative contributions within the expo- 
nential overlap. This same simple reasoning should play 
a significant role in our intuition about effective opera- 
tors of even higher mass dimension which are constructed 
from light fermion fields. The more light fermion fields 
that are present in the effective operator the more nega- 
tive contributions we can expect within the resultant ex- 
ponential wavefunction overlap. The significance of this 
result is that it shows that the geometric suppression 
from the warped RS background is sufficient to suppress 
the n-n transition operators to within the current exper- 
imental limits without any fine tuning while the effective 
mass scale can be as low as a fraction of a TeV. Any 
baryon number violating physics that may take place at a 
scale much higher than the electroweak symmetry break- 
ing scale in the 5D theory should therefore get warped 
down to the TeV scale in the 4D effective theory. 
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